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We investigate tlie evolution of gravitational wave perturbations about a brane cosmology em- 
bedded in a five-dimensional anti-de Sitter bulk. During slow-roll inflation in a Randall-Sundrum 
brane-world, the zero mode of the 5-dimensional graviton is generated, while the massive modes 
remain in their vacuum state. When the zero mode re-enters the Hubble radius during radiation 
domination, massive modes are generated. We show that modes decouple in the low-energy /near- 
brane limit and develop perturbative techniques to calculate the mode-mixing at finite energy. 



I. INTRODUCTION 

The Randall-Sundrum (RS) scenario [1, 2] is a sim- 
ple and novel way of realizing the idea that our ob- 
servable universe could be a 4-brane embedded in a 
higher-dimensional bulk spacetime, with interactions of 
the standard model confined to the brane, while grav- 
itational interactions access the bulk. The original RS 
model considered a Minkowski brane in an anti de Sitter 
(AdS) bulk, and this was generalized to a Fricdmann- 
Robertson- Walker (FRW) brane in a Schwarzschild-AdS 
bulk [3, 4]. 

The massless 5-dimensional graviton has a massless 
(zero) mode when projected onto the brane, but it also 
has a tower of massive Kaluza-Klein (KK) modes, i.e. 
modes which have an effective mass from a brane view- 
point. These massive modes introduce new features into 
the generation and evolution of cosmological gravita- 
tional waves, and the observational implications of these 
features provide in principle constraints on the RS sce- 
nario. It has been shown [6] (see also Refs. [8-10]) that 
during high-energy de Sitter inflation on the brane (in 
an AdS bulk), the zero mode is generated with a scale- 
invariant spectrum, at an amplitude which can be much 
greater than the corresponding general relativity value: 
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Here the general relativity expression is in square brack- 
ets, and the braneworld modification is in curly brack- 
ets, with £ being the curvature scale of the AdS bulk. 
(Table-top experiments to test deviations from Newton's 
law currently impose the constraint £ < OA mm.) At high 
energies, £H ^ 1, the correction factor is large {-^ ^£H). 

The massive modes are not excited during inflation [6] , 
and the zero mode remains frozen after inflation while it 
is beyond the Hubble radius. However, when the zero 
mode enters the Hubble radius during radiation or mat- 



ter domination, the separation between the massless zero 
mode and massive bulk modes no longer holds. A chang- 
ing Hubble parameter induces mode-mixing, and an ini- 
tial zero mode will generate massive modes. An estimate 
of this effect [8], based on an instantaneous transition 
from a de Sitter to a Minkowski brane, indicates that 
the effect will be very small. This estimate has been re- 
fined by considering the sharp transition to be from one 
de Sitter brane to another with slightly smaller Hubble 
rate [10]. 

In this paper we investigate the evolution of tensor 
metric perturbations about an AdS bulk in Gaussian 
normal coordinates defined with respect to an FRW 
brane. We first review the analytic results for perturba- 
tions about the original RS solution with two Minkowski 
branes embedded in AdS. We go on to develop a per- 
turbative analytic approximation (valid at low energies, 
near the brane) to study the generation of massive modes 
with an FRW cosmology on the brane. 



II. GRAVITATIONAL WAVE EQUATION AND 
ITS EXACT SOLUTIONS 

The generic 5D bulk metric that allows a spatially flat 
FRW metric on the brane, can be written in Gaussian 
normal coordinates as 

'^^W = -N^{t, y)dt^ + A^{t, y)dx^ + dy^ , (2) 

where the brane is located at y = 0. Tensor metric per- 
turbations [5-10] are given by 

(5)ds2 = ~N^{t, y)dt^ + A^it, y) [5y + h,j] dx'dx^ + dy^ , 

(3) 

where hij is 3-transverse (d^hij = 0) and 3-tracefree 
{S^-' hij = 0). We will treat one Fourier mode at a time, 
so that 
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with Cij a transverse and tracefree polarisation ten- 
sor which is an cigcnfunction of the spatial Laplacian 
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The tensor amplitude F obeys the 5D wave equation [6] 
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In the RS scenario, the bulk is an orbifold with Z2- 
symmetry about two fixed points, y = and y — L. 
This identifies points y ^ —y and y + L ^ L — y. We 
will assume that we have branes at the two orbifold fixed 
points with just the right surface energy-momentum ten- 
sors to satisfy the Israel junction conditions. In the case 
of Minkowski branes, this simply implies that the branes 
have constant brane tension. For a generic FRW brane 
at 2/ = this puts a strong physical restriction [12, 15] on 
the required equation of state on the second FRW brane 
at y = L, which wc will discuss in Sec. Ill A. 

Perturbations must also satisfy the Z2-symmetry. For 
a Z2-symmetric background this will be true, so long as 
the initial conditions are Z2-symmetric. The boundary 
conditions for the metric perturbations at the branes, in 
the absence of anisotropic pressure perturbations on the 
branes (which we will assume) are 



F'\y = O = 0, 



F'\y=L=0. 



(6) 



This should be imposed on any initial conditions, and is 
then preserved by the subsequent evolution. 

In the special case when the FRW metric reduces to 
a static Minkowski brane, the anti-de Sitter bulk metric 
Eq. (2) is given by [1] 



N = A = ex.p{-ny) , = 7 , 



(7) 



where /i is the mass scale corresponding to the anti-de 
Sitter curvature scale £. Here and in the following, the 
expressions hold for < y < L, and we consider the 
Z2-symmetry as implicitly imposed for other values of y. 
The wave equation (5) reduces to the simple separable 
form, 



where c±,ci,C2 are constants and B2 is a linear combi- 
nation of Bessel functions of order two. We only want 
the KK modes that satisfy the boundary conditions (6). 
For the zero mode, this requires 



C2 = 0. 
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For the massive modes, the first condition (at y = 0) 
requires 
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where the prime here is the derivative^ with respect to 
the argument of the Bessel function. This chooses the 
appropriate linear combination of J2 and Y2. The second 
(at y = L) requires 
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These two conditions can be rewritten, using the recur- 
rence relation for Bessel functions xB'^ + nBn = xBn-i, 
in the simpler form 



(18) 



where Bi is the same linear combination as B2 but with 
Ji and Yi instead of J2 and ¥2- 

The second condition can only be satisfied by particu- 
lar values of m^, which selects a discrete set of KK eigen- 
modes in the two-brane scenario. The allowed values of 
m are the solutions of 



Ji 



, m 
Yi - 



0. 



e'^'^y (f + k'^F^ = F" - AnF' . (8) 
Separating variables as 

F{t,y) = Y,Vm{t)Fm{y) , (9) 

m 

we find that 

(Pm + irri^ + k'^)^m = 0, (10) 
K-4/.F;, + mV''^i^ = 0. (11) 

The general solutions are [1] 

= c+e+*^"'+'='* + c_e-^^'"'+'='* , (12) 
Fo = Cl+C2e4''^ (13) ^lie'^e 

Fm = e^'^^Ba (^^ e'^y^ (m > 0) , (14) 



(19) 

The only case other than Minkowski branes for which 
the 5D wave equation (5) is separable in Gaussian normal 
coordinates is the case of a (anti-)de Sitter brane, which 
reduces to the Minkowski brane as a limiting case. For 
completeness we give the de Sitter mode functions in an 
appendix. 



III. COSMOLOGICAL BACKGROUND 
SOLUTION 



The 5D bulk metric, Eq. (2), for a general FRW brane 
at 2/ = in an anti-de Sitter bulk is given explicitly by [3] 

A{t, y) = a{t) {cosh /it/ -[1 + r{t)] sinh/it/} , (20) 

N{t, y) = + + wit)] r{t) sinh y.y , (21) 
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Here A is the tension for a Minkowski brane, so that the 
high-energy regime is p > A, and = Sn/M'^, where M5 
is the fundamental Planck sc ale. T he effective Planck 
scale on the brane is Mp = sj M| l\x. 

The two unknown functions of time, aif) and r(f), are 
obtained by solving the Priedmann and energy conserva- 
tion equations on the brane, 



a 



V/"^r(2 + r) , 
r = -3il(l + «;)r. 



(23) 
(24) 



For r <C 1, one recovers the standard Friedmann equa- 
tion. For r > 1, the evolution is unconventional. 
For a barotropic fluid with constant w we have [3] 
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(25) 

3(1 -h «;)/ui,(26) 



where a = a* when r = 1 , so that a* marks the transition 
between high- and low-energy regimes. In the case of 
radiation domination, w 



|, we have 



a = a, (l + 2/ii)^/^(4^f)i/S 



(27) 
(28) 



recovering the standard evolution, a oc t^^^, for t ^ A*~^- 
The coordinate singularity Pc, defined by A{t, Uc) = 0, 
and the event horizon yh, defined by N{t,yh) = 0, are 

given by 

coth /ij/c = 1 + r , coth /iy/i = 1 — (2 + Sw) r . (29) 



|, whereas yc exists 
00 as r 



There is no solution for y^ ii w > 

for all p > branes, and yc ^ 00 as r — > 0. For — 1 < 
w < —2/3 wc have yc < Vh, while for a de Sitter brane 
{w = -1), we have yc = yh- 



A. Regulator brane 

To simpliiy the problem, and make it numerically 
tractable, we will assume that there is a second regu- 
lator brane at a fixed position y2 = L < yc- Since ydt) 
increases with the expansion of the universe, the second 
brane remains within the regular coordinate region. The 
horizon y^ only exists for w < — | , which is the condition 
for inflation in the high-energy limit [11]. 

The flxed position of the regulator brane constrains its 
energy density and pressure via [12] 



Pl + A2 sinh /zl/ — (1 + ?") cosh fiL 

A cosh fiL — {1 +r) sinh fxL ' 

PL - A2 _ 2{pL + A2) 
A ~ 3A 
1 f sinh/i_L — [1 — (2 + 3iv)r] cosh fiL 
3 \ cosh/uL — [1 — (2 -|- 3u')r] sinh/xL 



(30) 



,(31) 



These imply 

(1 -|- WL)rL = —(1 + w)r [cosh/ui — {1 + r) sinh/xL]~^ 
X {cosh nL-[l-{2 + 2w)r] sinh M^}"^ , (32) 

where = p^/X, = p^/ p^. It follows that if the 
physical brane is de Sitter {w = —1), then so is the reg- 
ulator brane (with wl = —1)- For w > — | we have the 
limiting cases. 



r > 1 : {I + wlVl- 



[l + w) 



(33) 



[(2 Sw) sinh^ iiL\r 
r <Cl : (1 + wlYl « (1 + w)re^^ . (34) 

Thus the regulator brane requires (1 + wl)'I'l ^ for 
both the very high energy regime on the physical brane 
r ^ 00, and the low energy regime on the physical brane 
r ^ 0. 



B. New coordinates 

It is convenient to introduce a new dimensionless vari- 
able for the normal coordinate, 



z = ei^y , l<z<Zc 



1 + ^ 

r 



(35) 



where the range of validity follows from Eq. (29). The 
metric then reads 



-N'^it, z)df + A^{t,z)dx'^ + 



with Eqs. (20) and (21) being rewritten as 



A{t,z) = 



a{t) 



(36) 



(37) 



N{t,z) = !^|l+[2 + 3^i)]^[-'-l]}-(38) 



The wave equation (5) becomes 
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(39) 



where the prime from now on denotes d/dz. 



C. Causal propagation in the bulk 

The regulator brane will eventually have an effect on 
the physical brane via gravitational wave propagation. 
Approximate and numerical results will only be reliable 
up to the time when perturbations from the regulator 
brane reach the physical brane. The fastest signals in 
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the bulk are zero-momentum (k = 0) on the brane. They 
follow null world-lines with (taking n = 1) 



dz 
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(40) 



where, by Eq. (35), r{z^ - 1) < 2. Hence for a FRW 



brane with w > — |, we have 



< 3{l + w)ti. 



(41) 



In particular there is a lower limit on the time for causal 
propagation from the regulator brane at 2; = 2;^ = e^^^ 
to our FRW brane at 2 = 1, given by 



3At(l + w) 



(42) 



IV. 



COSMOLOGICAL BRANE 
PERTURBATIONS 



In general, the wave equation (5) for an FRW brane 
does not separate in a Gaussian normal coordinate sys- 
tem where the physical brane is fixed at y = 0. The 
wave equation is only separable if the background metric 
functions A{t, z) and N{t, z) arc themselves, and it can 
be shown that this is only possible in the special case of 
constant brane surface density (which is the case for a 
Minkowski brane [2] or de Sitter brancs [6, 13]). 

Moreover, even if one could define a set of independent 
eigenmodes for the bulk spacetime, the regulator brane in 
general represents a time-dependent boundary condition. 
Thus the regulator brane, moving relative to the normal 
coordinate system of the physical brane, will also lead to 
mixing between modes. In what follows we assume that 
the regulator brane remains at a fixed normal distance to 
the physical brane. One would expect this to underesti- 
mate the mode-mixing obtained when the second brane 
is free to move. 

The form of the bulk metric functions given in 
Eqs. (37) and (38) suggests that it should be possible 
to obtain an approximate separable solution in the low- 
energy limit, r — > 0, or close to the brane, 2; — > 1, where 
A and N take the limiting forms 



A(M)-^, 

Nit.z)^^. 
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In this limit we can study the limiting behaviour of F an- 
alytically. (A similar approximation was independently 
discussed in Ref. [14].) 

Formally we define the low-energy /near-brane regime 
by the condition 



r(22 - 1) < 1 . 



(45) 



If the regulator brane is fixed at finite z^ < Zc, then at 
sufficiently late times (in an expanding cosmology with 
w > —1) we will have r <C i/{z^ — 1) and the asymptotic 
solution will become a good approximation throughout 
the (finite) bulk. Conversely, the low-energy /near-brane 
condition, Eq. (45), will inevitably breakdown near the 
regulator brane in the limit Zl — > Zc- 

We can analytically estimate the effect of mode-mixing 
at finite (but still small) r by trying to build up a pcrtur- 
bative solution starting from the low-energy/near- brane 
solutions, presented in the next section, and then calcu- 
lating the corrections 0{r), which we go on to do after- 
wards. 



A. SepEirable solution for low-energy/near-brane 
limit 



Using the canonical variable, 

F = az-^l'^F , 



(46) 



and conformal time t] (so that n = a), we can write the 
wave equation (39) for r{z'^ — 1) — > as 

J_p + ^#(0) - ^2p^i?(o) = , (47) 
where the second-order sclf-adjoint operators are 

drf a ' 
^2 15 



9^2 4^2 ■ 



(48) 

(49) 



Here and from now on, a dot denotes d/drj. 

The wave equation (47) is separable and hence its so- 
lution can be expressed in terms of KK eigenmodes. 



F^^\'n,z) = Y,v'^^\ri)i,rn{z), 



(50) 



where the mode fimctions are related to those in Eq. (9) 

by v'^m = aiprn and i/'m = J^mZ~^^'^Fm, with A/to a nor- 
malization constant. They obey the equations 

(52) 



m 



'Dzfpm = o'^m ■ 

The bulk eigenmodes obey the same equation as in the 
case of a Minkowski brane, Eq. (11), and hence we have 
(for m > 0) 



V'm = J^m\/zB2 ( — Z 



(53) 



where the coefiicient of normalization ensures that the 
V'm constitute an orthonormal basis, 

dz tpm{z)tlJm'iz) = Smm' ■ (54) 



r 
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In the radiation era, w = ^ , Eq. (28) shows that in the 
low-energy/ late-time regime, 



and Eq. (51) becomes 



(55) 



(56) 



On large scales, or at late times, we can neglect the k- 
term, and the solutions are: 



.,(0) _ 



= CiT) + C2 



,(0) 



(57) 

(58) 



It follows that the massive modes m decay on super- 
horizon scales, unlike the massless mode m = 0. For 
T]^ 1 and k negligible, the massive modes behave as 



„(0) 



-1/2 
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B. Perturbative mode-mixing 

Using the bulk solution given in Eqs. (37) and (38), we 
can now write the wave equation (39) to first order in r 
as 



(60) 



S[F] = 



(2 -I- 3w)Vr,F 



-{l+w){5 + Sc1)n{F- - UF) - k'^F 



+ (1 - iw) 



zF' + If 



(61) 



where H = a/a = dj^aja and = p/ p. The full solution 
is a series expansion. 



p = ^(0) ^ ^(1) ^ _ _ _ 



(62) 



where the zero-order solution F^^^ is given by the solu- 
tion of Eq. (47) and successive terms in F correspond 
to successively higher-order terms in r in the wave equa- 
tion (39). Higher-order corrections can themselves be 
given as a sum over the zero-order bulk eigenmodes: 



F«(r?,^) = ^z;W(,,)V,„(^). 



(63) 



Being eigenmodes of the self-adjoint operator Dz in 
Eq. (49), the bulk modes ipm form an orthonormal basis 
for any function. 

The first-order corrections are given by 

Ip^/'d) + ^F(i) _ ^2p^/,(i) ^ rS[F(o\ (64) 



Substituting the decomposition (63) and projecting onto 
the basis of the functions leads to the following equa- 
tion that each of the coefficients v'^ {q) must satisfy: 

= r J ^ (/„„ - 5mn) [-(2 + 3w)n\i0) 

I n 

- A(i + ^)(5 + 3aH{^(°) 



2a2 

-3(l+«;)^^;^ 



(0) 



+ (l-3w;)5^( J, 



3. 



(65) 



Here the matrix coefficients and Jmn are given re- 
spectively by 



(66) 



Jmn = J dzipm{z)z — ilJn{z) . (67) 

It is the presence of off-diagonal terms in the matrices I 

and J that lead to mode mixing for r > 0. Note that for 
n = we have ipo oc z~^^'^ and hence Jmo = —3/2SmO- 

For a regulator brane at fixed zl, we expect mode- 
mixing to occur mainly at early times (maximum r) and 
become small for r <C l/(z|, — !)• Both analytic and 
numerical approaches are thus limited to finite < Zc- 
At a given time (fixed r) the mixing becomes largest at 
large z. However the coordinate si ngularity necessarily 
limits our analysis to z^ < Zc = -^/l + 2/r, and hence 
r{zj^ — 1) < 2 in all cases, suggesting a perturbative anal- 
ysis should not be too bad in most cases. 

The zero-mode growing-mode solution for k = given 

in Eq. (57), has vj^^ = Hvl^^ and hence the whole right- 
hand-side of Eq. (65) vanishes even at finite r. However 
an initial zero-mode configuration (v„i = for all m ^ 0) 
with finite k evolves into a mixed mode solution because 
the zero-order Vq^^ ^ acts as a source term at first order 

for all Vm^ with Imo 7^ 0. 

We can see from Eq. (65) that the zero-mode with 
n = evolves independently of the massive modes only 
when r = (Minkowski brane) or w = — 1 (de Sitter). 

We have checked the approximation by numerical so- 
lution of the full wave equation (39), in the case where 
w = ^ and only the lowest eigenmode (m = 0) solution is 
excited at lowest order, corresponding to an initial form 
of F that is constant in z, and with F = 0. In this case 
all the massive modes vanish initially. 



V. DISCUSSION 

We have shown how in the low-energy/ near-brane 
limit we can decompose bulk metric perturbations into 
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RS modes which evolve independently in the low-energy 
limit (r — > 0). At finite r > and away from the 
brane, z"^ > 1, the bulk metric is not separable in Gaus- 
sian normal coordinates defined with respect to a generic 
FRW brane and this leads to mode mixing. At late 
times/low energies, the bulk metric becomes separable 
for r{z^ - 1) < 1. 

Our analytic approximation shows explicitly how 
mode-mixing occurs when an initial massless mode, gen- 
erated during inflation, re-enters the Hubble horizon. 
The key result is Eq. (65). 

On the physical brane at z = 1, the 4D tensor metric 
perturbations, which are in principle constrained, e.g., 
by cosmic microwave background observations, have an 
amplitude 



1 



brane 



a(ry) 



fW(,?,1) + F«(,7,1) 



(68) 



to lowest order in our approximation. The massive modes 
contribute to the 4D tensor metric perturbations but 
their amplitude is suppressed at the brane due to the 
RS volcano- type potential for < < 15/^^/4 [2]. The 
massive modes contribute an anisotropic stress term in 
the "dark radiation" term [7] 



~A 



(69) 



where the derivatives here are with respect to y. 

We note that any analytic or numerical analysis based 
on Gaussian normal coordinates faces some significant 
practical limitations that may mean that it is not par- 
ticularly well-suited for more detailed calculations. In 
particular, the existence of a coordinate singularity at 
z = Zc, a finite proper distance from the brane, makes it 
impossible to treat an infinite bulk. Introducing a second 
brane at a fixed Gaussian normal distance z = zl < Zc 
may require an unphysical equation of state on the sec- 
ond brane. Thus unphysical effects may propagate to the 
"physical" (e.g., radiation-dominated) brane in a finite 
time. 
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APPENDIX A: DE SITTER BRANE 



The bulk background metric is given by 

(5)ds2 = N^{y) [-df + e^""'5^,dx'dx'] + dy^ , (Al) 
where a = e^"*, Hq is constant and N{y) is given by 

N{y) = cosh fiy — {1 +r) sinh /xy , (A2) 
with r constant. The wave equation has the form 

F + 3HoF + ^F = N'^F" + 4NN'F' , (A3) 



7 



and is separable, so that F{t,y) can be written in the 
form of Eq. (9), and 



and 



ifim + ^Hoipm +{m'^ + ^]ipm = , (A4) 



. (A5) 



N' rn^ 
N N"^ 

The general solutions are given by 



Fm{y) 



A. 



3/2 



7V3/2 



(A8) 



(PmW = exp ( --^Hot ) 



-Hot 



Ho 



with 



2 y m 



(A6) 



(A7) 



where ^13^2 ^ linear combination of associated Legendre 
functions. Equation (A7) indicates the existence of a 
mass gap [6, 9, 13] between the zero mode and the start 
of the massive KK tower at m ^ 



